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A new mechanism of particle acceleration, driven by the rotational slow down of the Crab pulsar, is explored. 
The rotation, through the time dependent centrifugal force, can efficiently excite unstable Langmuir waves 
in the electron-positron (hereafter e*) plasma of the star magnetosphere. These waves, then. Landau damp 
on electrons accelerating them in the process. The net transfer of energy is optimal when the wave growth 
and the Landau damping times are comparable and are both very short compared to the star rotation time. 
We show, by detailed calculations, that these are precisely the conditions for the parameters of the Crab 
pulsar. This highly efficient route for energy transfer allows the electrons in the primary beam to be 
catapulted to multiple TeV (~ 100 TeV) and even PeV energy domain. It is expected that the proposed 
mechanism may, unravel the puzzle of the origin of ultra high energy cosmic ray electrons. 

Conversion of even a small fraction of the enormous energy released in the rotational slowdown of pulsars 
(spinning neutron stars) into particle kinetic energy could create ultra high energy cosmic ray electrons. 
Most significant observational results in the realm of high energy cosmic ray electrons, beginning to be 
explored by the H.E.S.S., Pamela and Fermi collaborations' *, maybe summarized as: 1) the 2008 announcement 
of the H.E.S.S. collaboration showing evidence for a substantial steepening in the energy spectrum of cosmic ray 
electrons above 600 GeV, 2) Detection of cosmic-ray electrons up to 1 TeV, and confirmation of an excess of 
VHE leptons in cosmic-rays'^''' by the new instruments on board the Fermi spacecraft have. A thorough analysis of 
the observational data (of the Fermi telescope) on millisecond pulsars (the surface magnetic field close to B^i 
~ 10" "" G) by Kisaka & Kawanaka points to the existence of extremely high energy (s 50 TeV) electrons^. 

The pulsar rotation, and the dark matter annihilation have been suggested as two possible energy sources that 
could accelerate electrons to such ultra-relativistic regimes'*. All known pulsars are observed with decreasing 
spinning rates. The Crab pulsar, rotating with a period P ~ 0.0332 s, has a characteristic slow down rate 
P = dP/dt~4.21 X 10^^'ss^^ The corresponding slowdown luminosity, W~IQdQ/dt, is determined by the 
moment of inertia I = 2MR^^/ 5, the star radius R^f —10'' cm, the angular velocity Q = 2n/P, and the pulsar mass 
M~ 1.5 X Mq where Mq =2 x lO^'g is the solar mass. Out of this prodigious output W=5.5 x 10'*erg/s, only a 
part ~ 2 X 10^'" erg/s is sufficient to power the Crab Nebula; the rest - almost 60% of the released energy is, thus, 
available and seeks new channels of exploitation. 

The substrate for affecting the desired energy transformation is provided by the star magnetosphere. In the 
standard model for the pulsar magnetosphere"'^, an electrostatic field uproots electrons from the pulsar's surface, 
and accelerates them to relativistic energies along the curved magnetic field. The relativistic electrons emit 
copious curvature radiation. The emitted photons that have energies beyond the pair production threshold 
2 mc^ (m = electron mass ,and c = speed of light), in turn, create electron-positron pairs via y + B ^ + 
+ y' (B is the dipolar magnetic field). Newly produced pairs are further accelerated and emit photons; The 
cascade lasts till the nascent pair plasma, especially in regions farthest from the star, is dense enough to screen out 
the initial electrostatic field"''. 

The height h of the gap region (where the electric field is nonzero) is, therefore, strictly limited imposing an 
upper limit on the maximum attainable electron energy. The potential difference in the gap has been estimated to 
be A V=»l .6 X 10'^Bi2'^^P"^/^V'°, where B12 = -B/10'^ G ~ 6.7 is the dimensionless magnetic induction. Thus in 
the potential difference predicted by the standard gap model, the electrons may reach Lorentz factors of ~ 4 X 1 0" 
- not sufficient to explain the observed very high energy (VHE) radiation. Figure 1 shows the schematic repres- 
entation of the distribution function of plasma particles in standard models of the pulsar magnetosphere. The 
primary beam particles are shown by the narrow shaped area with the Lorentz factor yi,. The wider region 
represents the distribution of secondary electron-positron pairs, characterized by a corresponding Lorentz factor, 
yp when the distribution function has a maximum value. 
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Figure 1: | The distribution function versus the Lorentz factor. As is clear 
from the plot, the function consists of two parts: the first - wider region 
concerns the plasma component corresponding to the cascade processes of 
pair creation and the second one characterizes the primary Goldreich- 
Julian beam electrons. 

Results 

In this paper, we have constructed a theoretical framework for a very 
efficient and powerful mechanism for generating ultra high energy 
cosmic ray electrons. A highly relativistic electron-positron plasma 
in the magnetosphere of a pulsar, converts, the pulsar rotational slow 
down energy into Langmuir waves at rates that are several orders of 
magnitude faster than the star rotation period. These waves, in turn, 
via the well known and understood process of Landau damping, 
rapidly and preferentially, transfer the electric field energy to very 
high energy primary electron beam, accelerating them to even higher 
energies. It is shown that the radiation losses of various kinds do not 
much affect the very fast energy transfer process. Since the pulsar 
slowdown energy budget is enormous, charged particles could draw 



almost unlimited amount of energy via this mechanism. This rather 
general mechanism, that could encompass, for instance, the magne- 
tospheres of active galactic nuclei (AGN), and accelerate charged 
leptons as well as hadrons, could crucially advance the understanding 
of cosmic ray energetics. 

Discussion 

The determining parameters of the system are: the ratio of the mode 
frequency to the star rotation frequency w = (a JO., and 
0! = co2^/cJi = y-^ ly\^ and h. From our extensive numerical studies, 
we pick up two typical Mathematica plots (Fig. 2 and Fig. 3) in which 
ReNi, ImNj, ReNj, and ImN2 are plotted as functions of time. Time 
has been normalized to the inverse of the frequency cOi, the relat- 
ivistic plasma frequency corresponding to the stream with lower y 
(and much higher ambient density). The parameters for the illustrat- 
ive cases are: For Fig. 2,w = b= 10', a = .01, and for Fig. 3,w = b = 
10", a = .01. Since the x axis denotes time measured in their respect- 
ive plasma time (inverse of cUi) - the "absolute" time scales in the two 
figures are different by a factor \/lO. For each case, we choose initial 
conditions (t = 0) for what we call the Sin solution: Nj (0) = 0, ' (0) 
= 1,]V2(0) = 0,]V2'(0) = O!(cos -! sin b); it is easy to verify that these 
constitute a consistent solution to the coupled differential equations 
in the limit t^O. 

Here are the principal relevant results of our numerical investi- 
gations: 

a) After an initial gestation time t/w 1 (not displayed in the 
figures), the instability builds up with substantial growth rates 
as b approaches w; the latter is the regime in which the two 
streams interact, the interaction being the cause of instability. 
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{Sin, lr=S, li=Dot, 2r=BDash, 2i=SDash} 
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Figure 3: | The temporal behaviour of amplitudes of ReJVi (S), ImN^ (Dot), ReN2 (BDash), ImN2 (SDash). 



b) Since the two streams are assumed to have different intrinsic 
plasma frequencies (!X = cOj/cOj =yi'*/y2^1) the time 
dependence of (ReNi, ImNi), and {ReN2, ImNj) are quite 
different. It can be seen from the figures that although 
(ReNi, ImNi) oscillate (at a new slow frequency, the plasma 
oscillations are very fast and ever present but are not resolved 
in the pictures) and grow, the amplitudes {ReN2, ImN^) sim- 
ply grow. The growth is exponential as predicted by the 
Floquet theory of differential equations with periodic coeffi- 
cients. For the solution in Fig. 2, the growth time F"' ~ 2000 
plasma times = 2 X lO""* P, while the growth time for the 
solution in Fig. 3 is F"' = 6000 plasma times = 6 X 10"^ P. 
The growth rates (times) are very high (small) compared to 
the frequency (period) of star rotation. 

c) For a moderate and not too small a, the primary determinant of 
the instability gestation time, and growth rate is the relative 
value of b = 2 cklD. sin and w = oji/fl. Out of these, cui 
and are determined by the pulsar phenomenology. To bring b 
to the range of strong instability, one has to choose fc sin 0_ to be 
in the desirable range. In principle, then, for any finite sin 
one can find a range of k for strong growth. The centrifugal drive 
for Langmuir waves is, therefore, very robust and effective over a 
wide spectrum. 

The next major step, after establishing that the star rotation energy 
can be rapidly converted into charge separation electrostatic fields, is 
to show how this field energy can be harnessed, via Landau damping, 
to boost highly energetic particles to super high energies. We remind 
the reader that, although, the Langmuir waves are produced and 
sustained by the bulk e* plasma characterized by the relatively lower 
y region of the distribution (Fig. 1), we want them to damp on the 



primary beam electrons that have Lorentz factors as high as }■;, — 10^ 
(small blip on the high y of Fig. 1)". 

Electrostatic waves in relativistic e* plasmas'"''" allow a very sim- 
ple solution when the wave phase speed is very close to the speed of 
light, Vph = cOrIk ~ c. In this limit, of utmost relevance to the energy 
transfer mechanism, the real part of the dispersion relation reads'"" 

co = kc-P{k-kQ), (1) 

where P = (>')/(2y') is a very small parameter () denotes momentum 
averaging, and kg is the limiting value of k for which the phase speed 
is exactly c. The condition k s ko, implying Opi, £ c defines the regime 
in which such a wave could strongly interact with physical particles. 

Since the Langmuir waves with k approaching kg, have phase 
velocities asymptotically close to c, they will preferentially damp 
on, and energize the primary electron beam particles endowed with 
the highest y. For optimal transfer of centrifugal energy to particle 
acceleration, the rate of generation of Langmuir waves (measured by 
F, the instability growth rate) and the Landau damping rate'"* 



«r/i 

should be comparable. In Eq.(2), n, y and ojp = a/ 4ne^n/m denote 
respectively the number density, the Lorentz factor and the plasma 
frequency of the primary beam. 

A set of representative two-stream unstable e* plasma parameters 
that yield F ~ F^o. is given by y2 = 4700, and if y^ ~ 800. The 
corresponding instability time= damping time of 10"* s is several 
orders of magnitude less than the pulsar rotation period; the rate of 
energy transfer to the waves, and then on to the particles is extremely 
fast. It is worth noting, that despite the extremely efficient energy 
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transfer to particles, the strength of the pumping process is main- 
tained because the Langmuir waves are continuously generated 
drawing energy from the star slow down; the latter is taken to be 
an "inexhaustible" source. The non-linear effects (caused by the 
interaction between waves and particles) should be investigated at 
some stage, but we do not expect substantial difference in the results 
due to the "constancy" of the generation process. 

Since it is the electron-positron plasma that supports and feeds the 
two-stream instability, the corresponding characteristic energy in the 
unstable waves is of the order of the energy accumulated (transferred 
from star rotation) in the plasma components. Remembering that 
the number density associated with the e*-plasma far exceeds that of 
the beam component, the energy gained per beam electron could be 
extremely large if a substantial fraction of total "available" energy in 
plasma components were transferred via the electrostatic fields. 
Dispensing with technical details, we state here that, at the end of 
the process, a beam electron could end up gaining energy as high as, 

E-^ . (3) 

where _F ~ 2 mcQ^(r)"^ is the centrifugal force. For typical mag- 
netospheric parameters (the stream density ~ 2.5 X 10" cm"^, 
the primary beam density n\, = Ugj ~ 2 X 10^ cm"^), the energy 
acquired by a beam electron could get to be very high - in 100 s of 
TeVs, even PeVs. 

Highly relativistic particles will, inevitably, loose energy via radi- 
ation. We have, therefore, investigated how radiation will affect the 
overall energy transfer process worked out in this paper. We find that 
neither the inverse Compton cooling (in Klein-Nishina regime) 
nor the synchrotron mechanisms impose significant constraints on 
the maximum available energy. Very high energy electrons, with 
particle energy density exceeding the magnetic energy density by 
many orders of magnitude, move, practically, in a straight line. For 
them, the curvature radiation is also not expected to interfere with 
our fast energy transfer mechanism. Equally important is the fact that 
the ultra-relativistic electrons will have extremely small interaction 
cross sections with pulsar wind particles and interstellar protons so 
that they can breeze through, almost unimpeded, through the super- 
nova remnants. 

Methods 

In this section we present our model for studying the energy pumping process from 
the rotator - pulsar. Here is a two-step theoretical pathway that can, via well- 
understood plasma mechanisms'^"^^, convert the pulsar slowdown energy to create 
super relativistic cosmic ray electrons. Taking the Crab pulsar {B^t ~ 7 X 10^^ G) to be 
a representative, we show that: 

1) The time dependent centrifugal force (exerted by the star rotation) can, para- 
metrically, drive electrostatic plasma (Langmuir) waves in the relativistic e~ 
plasma residing in the star magnetosphere. The Langmuir waves grow through 
a variant of the two-stream instability with growth rates several orders of 
magnitude faster than the star-rotation frequency. Unlike the so called oblique 
modes discussed in^", we focus on the electrostatic waves, because due to extre- 
mely strong magnetic fields in the pulsar magnetosphere the particle distri- 
bution is one dimensional. 

2) The electric fields carried by these efficiently excited Langmuir waves can, then, 
damp, equally efficiently (via Landau damping) on extremely energetic primary 
beam electrons to accelerate them to ultra-high energies. 

The bulk e~ plasma (relativistic) that sustains and amplifies the electrostatic waves, 
and the even more relativistic primary electron beam on which the waves damp, are 
both constituents of the same pulsar magnetospheric region. 

To develop an appropriate theory for the centrifugally induced parametric 
Langmuir instability, we assume the magnetic field lines in the magnetosphere to be 
almost rectilinear - a good approximation for distances less than the curvature radius 
of the field lines. In the rotating frame of straight field lines co-rotating with the 
pulsar, the metric is given by 



(4) 



d7 



E+ xB 



(5) 



/ — e ,p is the species index, e, is the charge of the corresponding particle, dx = i^dU 
C = (1 — n^r^/c^)^'^ is the so-called time lapse function,^= —V^jc, is the gravita- 
tional potential, = ( 1 — /c^ ) ^^'^ is the Lorentz factor, and p,- and V,- are, 
respectively, the dimensionless momentum (p; Pi/m) and velocity. The first term 
on the righthand side of the equation represents the centrifugal force, and plays a 
crucial role in our model. 

In the laboratory frame, Eq.(5) takes the appearance of a fluid equation (see 
method) 



dt 



+ (v,-V)p, 



= - c^y.cVf + -(E+iviXB), 



that, in conjunction with the continuity 

- +V-(n,Vi) = 0, 



dt 



and the Poisson equation, 



V-E = 47i^e,n,-, 



(6) 

(7) 
(8) 



defines the entire dynamical system. 

If the electromagnetic terms do not contribute to the equilibrium (throughout this 
paper equilibrium simply means the absence of electromagnetic fields) because of the 
frozen-in condition Eq -I- Voi X Bq — 0, Eq. (6) reduces to^^ 



0.\ 



2 /drV 



(9) 



from which we can calculate the particle trajectories. For highly relativistic particles {)' 
^1), the equilibrium equation of motion, controlled by the centrifugal force exerted 
by the star rotation, yields the simple periodic solution. 



r{t) = -^sm{ilt + (j>). 



Vo(f) = yosin(Qt + ii4), 



(10) 



(11) 

where Vq ~ c is the velocity amplitude. The phase <j), necessary to insure a general 
solution of Eq. (9), was taken to be zero in previous work^^. We wiU soon show that, 
due to this omission, a crucial ingredient of the current modela strongly growing 
variant of the "two stream instability" - was missed in the earlier work. 

In our model we are trying to describe an essentially kinetic plasma by a fluid 
approach. The particle distribution has a wide )' (energy) spectrum (Fig. 1), and time 
dependent particle velocities have a characteristic phase (j) (11). A proper detailed 
description is likely to be rather involved. In order to extract the main results in a 
highly reduced description, we will assume that the plasma is multi-stream; each 
stream has a characteristic y and a phase ^. By solving the linear, connected dynamics 
of two such streams, we will demonstrate the emergence of a very interesting mani- 
festation of the two stream instability; the instability will be determined by the phase 
difference between the streams. Each stream constitutes an electron-positron plasma, 
and Equation (6) is to be viewed as the equation of motion of the specie i (electron or 
positron) belonging to a given stream. 

The orbits (10,11) (calculated using the frozen-in condition), are time dependent 
and periodic with the periodicity determined by P ~ In/Q, the star rotation time. 
Strictly speaking, these are field-free but not equilibrium orbits; they may be con- 
sidered equilibrium orbits only as long as the characteristic wave time-scales assoc- 
iated with the wave motion (both oscillation and growth) are much shorter compared 
to P. Because of time dependence of the orbits, the linear wave equations will not be 
amenable to the standard time - Fourier analysis; the solution of the stability problem 
wiU be more involved and nuanced. The periodicity of the equilibrium orbits is what 
leads to the parametric pumping of particle energy into the energy in the electric field. 

Following the "standard" procedure of splitting the physical variables 4* — {n, v, p, 
E,) into equilibrium and perturbed parts, 4* ~ 4*° -I- expanding the linear per- 
turbations as (//^(f, r) cx )//^(0 exp [f(kr)], we derive, after some tedious but 
straightforward manipulations, the two coupled, non-autonomous differential 
equations 



d^N2 

"dp 



-(^2-^2^ 



in terms of the perturbed 'densities' of each stream 



(12) 



(13) 



(14) 



In the single particle equation of motion (in the I + 1 formalism^^) 



The streams 1 and 2 are distinguished through their characteristic y and the phase (j 
(for simplicity, we have put the velocity amplitude Vq — c). The frequencies 
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ftJi 2 = y Sne^ «i 2 / myl 2 respectively the effective relativistic plasma frequencies 
(for longitudinal modes) corresponding to the two streams, and 
ck 

— [sin(a^ + (/>i)-sin(at + ^2)]^^ cos (Of + (/> + ) , (15) 
with 2(j>± = (j)i ± 3-rid 

?7-2^sin <f)_. (16) 

In deriving Eqs.(12-13), we have invoked {ck/Q) ^ 1; this parameter wUl turn out to 
be approximately the ratio of the Langmuir frequency to the star rotation frequency 
and, will lie in the range 10^ - 10^". 

From the general structure of evolution Eqs.(12 - 13), with defmitions (15 - 16), 
one may draw the following inferences: 

1 ) The evolution of the two streams is coupled, and the time dependence of the 
coefficients is fuUy contained in 

2) Though / depends on both the sum and difference of the phases associated 
with the two streams, the dependence on the phase sum 2^+ = (!>i + <l>2 is 
trivial, and can be absorbed in a redefinition of the time variable; (j)+ will be 
dropped from now on. 

3) The phase difference <j)-, however, emerges as a fundamental parameter of the 
two-stream system. If (j)- ~ 0, then / ~ 0, and the coefficients of Eqs. (12-13) 
become constants (the system becomes autonomous) leading to the simple 
dispersion relation 0/ = co{^ + C02^, i.e, the two streams oscillate together with 
a frequency that is the sum of their individual plasma frequencies. Most 
importantly, for <j) = 0, there is no parametric instability. Strictly speaking, 
this is only an approximation, since in Eq. (15) we have assumed Vo/f ~ c, 

whereas Vo/(~c{ 1 ] , leading to a nonzero growth rate (even if ^_ — 0) 

which will be less than for non vanishing values of ^_ by many orders of 
magnitude. 

4) It is only when (j)- is non-zero (the two streams remain "distinct") that the 
system "notices" the periodic time dependence of the centrifugally driven 
orbits. Since the time dependence and periodicity of the orbits is, precisely, 
what drives the strong Langmuir-wave instability, assigning different phases 
to the particle orbits (10-11) is the key step that allows energy transfer (via 
instability) from rotational slowdown to the Electrostatic fields carried by the 
Langmuir waves. 

Differential equations with periodic coefficients (like (12-13)), Mathieu equation 
being the simplest example, are known to display parametric instabilities in well 
defined regimes. We find this behavior when (12-13) are solved, using analytical 
(asymptotic) as well as numerical techniques, to determine the time evolution of the 
densities Wi (and the electric field, £ — —iV2e'0- Although the analytic solution 
is quite instructive, we will dwell here, mostly, on the more general numerical solu- 
tion. Mathematica solutions demonstrate that model Eqs. (12-13), indeed, harbor 
parametrically unstable solutions. The charge separation densities (the electric field) 
grow in time and acquire a hefty growth rate when co/Q.; this qualitative feature is 
also predicted by the analytical dispersion relation. 
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